Electrostatic approximation of source-to-target mean first passage times on lattices 
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We demonstrate that the source to target mean first passage time (MFPT) is approximately given 
by the potential difference of an electrostatic problem which shows that the MFPT scales like the 
resistance between the target and a distant shell. This analogy allows the asymmetry of the MFPT 
on non uniform lattices to be incorporated and provides a number of useful insights. For example, 
on transient lattices, the MFPT converges exactly to the product of the mass and the site dependent 
resistance between the target and a shell at infinity. 



Problems which require calculation of the time it takes 
a random walker (RW) to reach one or more traps on a 
network arise naturally in numerous settings [l[ including 
chemical reaction rates 0, Q , transport in random me- 
dia [U, search algorithms, spreading of disease and other 
problems @ . The single most important characteristic 
of the density of trapping times is its average; the mean 
first passage time (MFPT). A variety of MFPT problems 
can be defined. These include the time it takes to reach 
a trap if the RW is released at a random site on the net- 
work [J H, H, H3] , the time to reach any site a distance r 
away (J| , or the time it takes a walker to visit a target x 
if it is released at a source y [1, [ll[ . The nature of the 
process to be modeled dictates the relevant MFPT. 

Recently, there has been interest in the scaling form 
for a variety of source to target MFPT problems on finite 
lattices p| [l2[ . Condamin et al derived the form 

( A + Br d ™- d f d w -d f >0 
(T(x\y))aN{ A + B\n(r) d w -d f = (1) 
[ A-Br d ™- d f d w -d f < 0. 

for lattices with N sites. Here d w and df are the random 
walk and fractal dimensions respectively. Interestingly 
A and B are independent of the lattice shape. The ap- 
proximation provides an extremely simple solution to a 
problem which one might expect to show a strong de- 
pendence on boundary shape and the source to target 
positions. 

Although Eq. (fTJ) was found to give very good predic- 
tions on a range of lattices, its derivation assumes that 
the diffusion propagator at all points of the lattice is iden- 
tical. This is only strictly true on "uniform" networks 
such as n dimensional lattices. Recent theoretical in- 
terest in the MFPT has been stimulated by the study of 
transport on complex networks [l3l . [l4| which have a high 
degree of non-uniformity in their coordination number 
and exhibit modularity. Moreover, there are numerous 
loopless networks which have non-uniform coordination 
numbers and branch lengths. The absence of loops is 
known to significantly affect their dynamics [l5j . 

It is therefore interesting to consider (r(x\y)} on non- 
uniform lattices in more detail. We do this by deriving 
a scaling form for the MFPT using an analogous elec- 
trostatic problem. Although our approach holds in gen- 
eral, we focus in particular on loopless networks and on 
classes of trees that do not follow the standard propa- 



gator form [16j |. The method is fruitful in the following 
respects; it incorporates the point to point character of 
the MFPT, leads to a more general scaling law, explains 
independence of the scaled MFPT ({r(x\y))/N) on lat- 
tice shape, and shows how t(x, y) is related to resistance. 

On a finite lattice it is possible to show that the MFPT 
r(x\y) for a walker released at y to arrive at x is exactly 
given by 0, [T3| 

r(x\y) = 2M [°° _ ^A dt (2) 

Here p(y\x,t) is the probability that a walker released 
at x at t = is at the point y at time t and M is the 
number of bonds within the structure. If the lattice is 
uniform (k x — K y — k) then 2M/k = N. The probabil- 
ity satisfies the discrete diffusion equation^ = V 2 p+<$t o 
with reflective boundary conditions which has steady 
state solution p(y\x,t) — K y /(2M). To obtain r(x\y) 
exactly, p has to be calculated for each finite lattice. In 
Ref. [5| the authors derive Eq. |T|) based on the hypoth- 
esis that p(x\x,t)/n x - p(y\x,t)/n y rj Poo{y\x, t)/n x - 
p oa (y\x,t)/ K y where p^ is the Green's function for the 
case where the lattice has infinite size. Certainly both 
expressions are identical for small times (before the re- 
bound from the boundary) and presumably have simi- 
lar long time behavior. Using a well known propagator 
form [H, [H[ p{x\y,t) w H(r dw /t)/t d f/ d ™ the integrals 
are evaluated directly, giving rise to Eq. ([TJ). As r(x\y) 
has a point to point character, the scaling law given in 
Eq. {!]) (which only depends on distance) refers to the 
MFPT to reach a target averaged over all sources a fixed 
distance away. We denote this averaging by (•). 

Instead of proceeding as above we show how the in- 
tegrals in Eq. ((2|) can be interpreted as electrostatic po- 
tentials. To see this consider terminating the integra- 
tion at a time scale associated with diffusion of the fi- 
nite lattice T = bL dw where I is a typical chemical 
length (i.e. stepping distance) on the structure and d w 
is the random walk dimension in chemical space. In 
Ref. [Tj| (also see Ref. (2C|) it is argued that the inte- 
gral P(y\x, T) = p(y\x,t)dt « K y (f>(y\x; L), where 4> is 
the potential at y due to the injection of a unit current 
at x with a Dirichlet boundary condition applied at L. 
P{y\x,T) exactly represent the concentration field cre- 
ated by the release of a walker at x at each of T + 1 time 
steps. The identification with the potential follows from 
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the hypothesis that P(y\x,T) has equilibrated near the 
origin (thus satisfying the potential equation), and that 
P(y\x, T) w on the boundary. The latter claim is plau- 
sible because only walkers released at small times will 
have reached d{x,y) = L, where d{x,y) is the chemical 
distance between x and y. 

Thus the MFPT is related to the potential by the ex- 
pression 

r(x\y) w 2MA0 = 2M[cf>(x\x] L) - </>(y\x; L)}. (3) 

To proceed recall that the point-to-shell resistance R x (£) 
is defined as the potential difference induced by the flow 
of unit current from the target x to a grounded shell 
at a chemical distance £. If the potential depends on 
the chemical distance £ alone then (A</>) = R x (£) which 
should be independent of L (the location of the Dirichlct 
condition). In general the potential depends on position 
y, but not too strongly for lattices which follow the fractal 
Einstein law [l6l |2l|. Thus, we find the approximation 

<T(x|y)> « 2MR X (£). (4) 

From the asymptotic form of the resistance for large 
I [13], the MFPT (for £ -> oo) will scale as 

( Bf C > 
{r{x\y)) w 2M I B\n(£) C = (5) 
[ R* - Bf C < 0. 

Here R* is the resistance from x to a shell at infinity and 
C is the resistance exponent. The resistance R x {£) is not 
expected to strongly depend on the size of the lattice as 
it is dominated by the shortest path(s) between x and y. 
Equation ([5]) therefore provides an important insight into 
two striking features of the MFPT identified in Ref. ||; 
the near linear dependence of the MFPT on mass, and 
the independence of the coefficients in Eq. ([T]) on the 
shape of the confining domain. 

If the fractal Einstein law £ = d w — dj is assumed, then 
the argument recovers the scaling form of Eq. ([1]) except 
that the coefficients are now determined by the electro- 
static problem. The electrostatic scaling forms (which 
are strictly true only for large I) give important insights 
into the constants A and B appearing in Eq. JT]). For 
( > 0, there is strictly no constant in the asymptotic 
form while for C < 0, NA = 2MR*. Note that there 
are lattices on which £ ^ d w — df, such as those created 
by diffusion limited aggregatio n 12311 (DLA) and certain 
asymmetric deterministic trees [la ]. Thus the resistance 
form expressed in terms of C is more general. 

The MFPT is also connected to the point-to-point re- 
sistivity p{x, y) which can be defined as the potential 
difference between two sites if a unit current is injected 
at one and withdrawn at the other. If the Laplacian 
of Eq. ([2]) is taken and the right hand side is integrated, 
r(x\y) can be shown to satisfy S/ y T(x\y) = l/M — 5 xy @. 
Now consider the equation satisfied by the function 
X(z) = [r(y\z) - t{x\z)]/{2M); V 2 zX = S xz - S yz with 



reflective boundary conditions. The latter equation is 
exactly the potential associated with the flow of unit cur- 
rent from x to y so pix, y) = x( x ) ~ x(v) which recovers 
a result of Chandra [24| 

T{x\y) + T(y\x)=2Mp(x,y). (6) 

Note that if r(x\y) = r(y\x) then the MFPT is exactly 
given by r(x\y) — Mp(x,y) (25| (see [26| for a point to 
shell analogue). This is obviously true on uniform net- 
works in which every site is identical (e.g. n-d lattices). 

In terms of current discussion, the fact that r(x\y) = 
Mp(x, y) on uniform lattices offers two interesting in- 
sights. First, we have argued that r(x\y) « 2MA<fi which 
implies that p(x, y) ~ 2A<fr. This result becomes exact for 
infinite uniform lattices and is in fact useful for calculat- 
ing p(x, y) [27[. Secondly, it allows a useful interpretation 
of the MFPT derivation given in Ref. [f|. For uniform 
finite lattices Eq.@ becomes p(x,y) = J °° p(x\x,t) — 
p(y\x,t)dt which should remain well defined for infinite 
lattices giving p M (i,y) = Poo(x\x,t) - Poo(y\x, t)dt. 
Thus the finite to infinite assumption can be viewed as 
r(x\y) = Mp(x,y) ~ A / Ip oc (x,y) which coincides with 
the form (g]). 

It is generally not clear how good an approximation 
r(x\y) — r(y\x) is for non-uniform lattices. Certainly it 
will hold exactly for pairs of points which are reflections 
in lines of symmetries of the structure (e.g. the Sierpinski 
gasket). An obvious counter example is if x is the end 
of a dead end branch and y is its neighbor. In this case 
r(x\y) = 1 whereas r(y\x) = 2M—1 since p(x, y) — 1. An 
expression for the asymmetric component of the MFPT 
is found by combining Eqs. ^ and ^ and using the 
identity n y p(x\y,t) = n x p(y\x,t) to obtain 

r(x\y) = Mp(x, y) + M f P J^A _ PjvM dt . (7) 

Note that the integral term of Eq. ([7]) has been studied 
in detail by Noh and Rieger [l7|. The formula is useful 
below, and also clearly demonstrates how the assumption 
of a site independent form of the propagator p(y\x, t)/ k v 
is connected to symmetric MFPTs. 

Benichou et al. [l2j give an exact expression for the 
average of r{x\y) over the k x neighbors of x as (r(l)) = 
2M/k x — 1. From Eq. (TJJ, this implies that the point 
to shell resistance R x (l) w (r(l)>/(2M) w l/n x . The 
approximation is actually exact because there 
resistors in parallel between the site x and the shell 
grounded at d(x,y) = 1. In Ref. [12 it was hypothe- 
sized that the constants A and B could be evaluated by 
using the boundary condition (t(1))/N rs 1 (these were 
termed 'zero constant' formulas [l2|]). If a similar analy- 
sis is performed for non-uniform lattices this would imply 
B = 1/k x for C > and B = R* - l/n x for ( < in 
Eq. (O . Note that the form of the resistance for large £ is 
asymptotic in nature and there is no guarantee that ex- 
trapolating the scaling form to £ = 1 to fix the constants 
will be accurate. We test the results below. 
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FIG. 1: The MFPT on a finite 3-d network (o) is almost 
exactly given by 2A(f> (+) on an "infinite" lattice 

If the MFPT is averaged over all sources and targets 
a chemical distance I apart then ?(£) = M p(£) exactly. 
Here p{£) is defined by the same averaging process (recall 
p is the point to point resistance). This can easily be 
seen from Eq. (J6J) or (JTJ) . For loopless networks p(£) = £ 
and therefore f(£) = Ml exactly. This form is confirmed 
by the numerical results in Fig. 2(b) of Ref. [l2j for the 
T-tree. However it is not clear that the coefficients of 
?(£) and (r(x\y)) are identical. Furthermore, the scaling 
exponents of averaged quantities may differ from that of 
the original quantity [28|. Indeed, for loopless network 
having £ < 0, the form {r(x\y)) ~ 2M(R* - Bf) does 
not scale as Ml. 

To investigate the approximation t ps 2MA<f> we have 
calculated the MFPT from each source on a range of 
finite lattices to a fixed target x. The MFPT r{x\y) 
is computed by releasing a walker at the target x and 
using the resulting p(y\x,t) to evaluate the integrals at 
each y of Eq. ((2]) exactly. The question is whether the 
potential difference Ac/> on an "infinite" lattice can give a 
good approximation to the MFPT on finite sub-domains 
around the target. We first consider the MFPT on a 
finite 3-d lattice I < 15 with a target at its center. The 
potential Acf> is computed on a 3-d lattice with Dirichlet 
boundary conditions at L — 50. In this case A(j> (Fig.[T]) 
reproduces the point to point MFPT almost exactly. 

The MFPT of a DLA cluster is shown in Fig. EH It is 
seen that A(f> captures many features of r(x\y). In par- 
ticular, the MFPT is seen to be nearly constant along 
dead end branches with value given by the MFPT at 
each branching point. This exactly matches the qualita- 
tive behavior of the potential. The solid line in the lower 
left box shows that (t) / (2M (Ac/))) « 1 (supporting the 
scaling form ([5])) but the point to point comparison be- 
gins to show significant outliers. The graph on the lower 
right shows that a practical approximation is obtained 
from t « M/M' x t' where the prime indicates a larger 
lattice. This is expected, and indicates that the MFPT 
on a lattices of size M < M' can be estimated accurately 
by multiplying r' by a simple factor. 

We now consider the three deterministic trees (DT- 
A,B,C) described in Fig. [3 The MFPT for DT-A (C = 
—0.58) is shown in Fig. 2] for the cases of two targets; 
one located at base of the tree and the second at its im- 
mediate neighbor. The MFPT exhibits mild anisotropy, 




FIG. 2: The MFPT of a DLA cluster clipped at £ = 80 (top 
left) compared with the potential on the same DLA cluster 
with Dirichlet conditions applied at £ — 140 (top right). The 
ratio of the two functions is shown bottom left. The final 
figure (bottom right) compares the scaled MFPTs obtained 
for the cluster clipped at radii £ — 80 & £ = 140. 




FIG. 3: An example of a deterministic tree (DT) constructed 
by attaching copies of the rescaled base unit to each end of 
the base unit's end points and so forth [l^.[29l]. Shown is DT- 
C which has (ui,li2) = (1,2). We also consider DT-A which 
has («i,«2) = (2, 1) and DT-B which has (u\,U2) = (1, 1) but 
the length of successive iterates quadruples. 

but is very well quantitatively approximated by (Ac/)) on 
a larger lattice. A similar picture is shown for tree DT-B 



which has £ 



The properties of both networks DT-A 



and DT-C obey the FE law ( = d w — df so follow the 
scaling of Eq. {l} (and therefore Eq. ([5])). In contrast, 
on network DT-C, £ ^ d w — df and transport is there- 
fore highly anisotropic. This is shown in Fig. [5l The 
point to point approximation A<p is only able to repro- 
duce the MFPT at small £ for this model, but it is clear- 
that (t) « 2M (A</>) (solid line, lower left box) confirming 




100 150 200 250 



100 „ 150 200 250 



FIG. 4: The MFPT for DT-A (left) (( = -0.58) for two 
different targets. The solid lines are (V(/>) and confirm that 
r(x\y) ~ 2MR X (£) with the limiting values given by R%. The 
results for the model DT-B (£ = 1/2) are shown on the right. 
For these models £ = d w — df. 
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FIG. 5: Results for model DT-C. The sub-plots match those of 
FigO For this lattice C, 7^ d w ~dj implying highly anisotropic 
transport. 
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FIG. 6: A comparison of (r) (o) to the zero-constant formula 
described in the text. In the left figure the dashed line is the 
asymptotically fitted form. The right figure shows the results 
for networks DT-B (□, solid line) and DT-A (o, dashed line). 

(t) w 2MR x and hence Eq. ©. 

In Fig. [5] we plot (r) at a single target vs the 'zero 
constant' formula with the values of B described previ- 
ously. On the uniform 3-d lattice the formula appears 



reasonable but if B is obtained from the large I behav- 
ior the value is 0.11 (dashed line in Fig. [6]) rather than 
0.0854. For the non-uniform DT-A and DT-B networks 
the approximations are seen to be poor and we conclude 
that the 'zero-constant' approach [12| does not generally 
work for non-uniform lattices. 

We have shown how the source to target MFPT on a 
finite lattice can be interpreted as a potential difference 
2AfA(/> on an infinite extension of the lattice. The ap- 
proximation 2MAtp reproduces the point to point charac- 
ter of r(x\y) of the function near the target, particularly 
for uniform lattices. There is some advantage in obtain- 
ing r{x\y) from a single calculation of A<f> in that an 
approximation for lattices of all sizes is obtained, how- 
ever a significantly better approximation is yielded by 
r(x\y) = (M I 'M')t> '(x\y) where r'(x\y) is the MFPT on 
a large lattice (bottom right panels of Figs. [2] and [5]). 

The spatially averaged MFPT (over sources at distance 
I) (t) 2M (A0) is identified as the point to shell resis- 
tance (t) 2MR x (£) because the potential <f> is associ- 
ated with a unit current flowing between the target and 
a distant shell. These expressions are confirmed by the 
numerical results. Note that the scaling forms ([5|) are 
more general than |T]) in the sense that they encompass 
lattices which do not follow the fractal Einstein law (e.g. 
DLA and DT-C). Moreover, the form does not involve the 
fractal dimension directly, and hence we propose that it 
can hold for lattices with a well defined exponent £ (e.g. 
bundled structures [30] such as combs.) 

Given the connection between the MFPT and resis- 
tance, we believe that further insights may be gained into 
other problems (e.g. a random walker leaving a confined 
domain jy]) by studying their electrostatic analogs. 
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